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Abstract 

We show that a passing gravitational wave may influence the spin entropy and spin negativity 
of a system of iV massive spin-1/2 particles, in a way that is characteristic of the radiation. We 
establish the specific conditions under which this effect may be nonzero. The change in spin entropy 
and negativity however, is extremely small. Here, we propose and show that this effect may be 
amplified through entanglement swapping. Relativistic quantum information theory may have a 
contribution towards the detection of gravitational wave. 



Relativity and quantum mechanics are the two pillars of 20th century physics. Einstein's 
general relativity |l| is the classical theory of gravity. Many of its predictions have been 
experimentally confirmed via very precise measurements. One of its most intriguing predic- 
tions is the propagation of ripples in spacetime curvature at the speed of light (c = 1) called 
gravitational waves. Gravity is a long-range interaction and it is not possible to shield this 
interaction. Gravitational waves thus provide a new window for exploring astronomical phe- 
nomena. However, gravity is the weakest of the four fundamental interactions; this means 
that gravitational waves are not easily detected. In fact, they have not yet been detected 
on Earth. Quantum mechanics, in combination with computation and information, leads 
to unexpected new ways that information can be processed and transmitted, extending the 
known capabilities in the field of classical information to previously unsuspected limits Q]. 
One of the greatest challenges faced by quantum information scientists is the fragility of 
quantum coherence and entanglement in the presence of environmental decoherence [3]. In 
particular, multipartite entangled states such as the GHZ states |4( become more suscepti- 
ble under certain kinds of noise as the number of particles increases |^|. Motivated by its 
fundamental importance in gravitational wave detection and several recent developments 
in relativistic quantum information theory p, we explore, in this paper, the possibility of 
turning this fragility into a quantum means to detect gravitational radiation. 

Peres, et al. |7] were the first to study the relativistic properties of spin entropy for a 
single, free particle of spin 1/2 and nonzero mass in flat spacetime. They showed that even if 
the initial quantum state of the particle is a direct product of a function of momentum and 
a function of spin, the state under a Lorentz boost is in general not a direct product. This is 
because the spin undergoes a Wigner rotation py whose direction and magnitude depend on 
the momentum of the particle. Spin and momentum appear to be "entangled". As a result, 
the reduced density matrix for spin becomes mixed and the corresponding entropy becomes 
nonzero. Slightly later, Gingrich and Adami jjj] showed that Lorentz boost can also affect 
the entanglement between spins. Namely, a maximally entangled Bell state of two massive 
spin-1/2 particles loses entanglement under a Lorentz boost. 

More recently, Terashima and Ueda |10J] (see also extended the original investigation 
by Peres et al, by considering the relativistic quantum mechanics of a massive spin-1/2 
particle moving in curved spacetime, which entails a breakdown of the global 5*0(3, 1) 
symmetry associated with flat spacetime. In this case, spin can only be defined locally 
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at each spacetime point by invoking the SO (3,1) symmetry of some local inertial frame. 
Specifically, a spin-1/2 particle in curved spacetime is defined as a particle whose one-particle 
states furnish the spin-1/2 representation of the local Lorentz transformation. Terashima and 
Ueda showed that, as a consequence of this local definition, the motion of the particle is 
accompanied by a continuous succession of local Lorentz transformations, which gives rise 
to spin entropy production that is unique to the curved spacetime. They illustrated their 
ideas with the Schwarzschild spacetime in jloL fll | . 

In this paper, we study the effects on the quantum states of one, two, and in general N 
spin-1/2 particles due to a plane gravitational wave spacetime propagating in the positive 
^-direction jlj|: 

ds 2 = -dt 2 + [1 + f{t - z)]dx 2 + [1 - f{t - z)}dy 2 + dz 2 . (1) 

The size and shape of the propagating ripple in curvature are determined by some dimen- 
sionless function / (\f(t — z)\ « 1). For example, for a Gaussian wave packet with width 
uj and maximum height A, 

(2) 

And, for a gravitational wave of amplitude A and definite frequency w, 

fit - z) = A sin[tz7(t - z)\. (3) 

The only nonvanishing Christoffel symbols for the above metric are 21^ = — 2Y yy = df/dt, 
-2Y% X = 2T* yy = df/dz, r* = dhi y/T+J/dt, T% z = dln^T+J/dz, Y y ty = d\n^T=J/dt, 
and Y\ z = 9 In y/T=J/dz. 

We begin with a single spin-1/2 particle A of mass m in a local inertial frame at the 
spacetime point xf . This particle is initially prepared at proper time in the state 

|V)a= / N(k a )d 3 kY,C(k a ,X)\k a ,\)A, (4) 

J A 

where Nik a )d 3 k = md^kj \Jk ■ k + m 2 is Lorentz-invariant volume element. From here on, it 
is assumed that Latin and Greek letters run over the four inertial-coordinate labels 0, 1, 2, 3 
and the four general-coordinate labels, respectively. \k a , A) a is the momentum eigenstate of 
the particle, labeled by the four-momentum k a = i\Jk-k + m 2 , k) and by the ^-component 



fit — z) = Aexp 



(t 
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A (=| or 0, I or 1) of the spin. We consider, in particular, the case where the coefficient 
C(k%\) = D(k a )5 X0 , 



D(k a ) 



n ex p 



(k a - q a {xi)f 
2^ 



(5) 



^/N(k a )^w o=i,3 

We assume that the spacetime curvature does not change drastically within the space- 
time scale of the wave packet. q a (xi) is as given in Eq.(18). Together with the or- 
thogonality condition A (k' a , \'\k a , \) A = S 3 (k' — k)5y\/N{k a ) we clearly have a(^\^)a = 
J N(k a )d 3 kJ2x\C(k a , A)| 2 = 1, i.e., \ip) A is normalized. To ease calculations, we set k 2 to 
zero with no loss of generality. It follows that at Tj the reduced density matrix for spin, 



PA ( n ) = J N(k a )d 3 k A (k a \^) A m a ) A = |0)a<0|, 



(6) 



and the corresponding entropy S A {ji) = — tr[p^(rj) log 2 pA( r i)] = 0. We will show that at a 
later proper time 77, p A (ri) evolves to 

\ 

(7) 



1 1 + c s 
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with spin entropy S' A {r f ) = -Plog 2 P - (1 - P) log 2 (l - P), P = (1 - \u\)/2. Here, 



u = J N{k a )d 3 k\D{k a )\ 2 exp{iQ), 



(8) 



with f2 = Q(k a ;Ti,Tf , £, 1?) as given in Eq.(22), c = Re(tt) and s = Im(-u). It will be useful 
to note 



£[iH EE £ 



£[R 10 ]=£[\l)(0\]= 1 - 



1 / 1+c s 

s 1 — c 

-s -1 + c 
1 + c s 



, £[i? 01 ]^[|0>(l|] 



-s 1 + c 
2 1 -1 + c s 



, £[ip] = £[ll)<ll] = - 



1 — c — s 
-s 1 + c 



• (9) 



Next, we consider two spin-1/2 particles A and B with equal mass m, initially prepared 
in the state 

\*)ab = J J ' N(k a )N(p b )d 3 kd 3 pY,C(k a A;p b ,v)\k a ,\)A®\p\a)B, (10) 

with C(k a ,\;p b ,a) = D{k a )D{p b )5\ a . By writing |\I/)ab as a density matrix and tracing 
over the momentum degrees of freedom, we obtain, at Tj, a maximally entangled Bell state 



XABin) = \* Bell) Ab(* Belli = ~ £ ^ ® 

Z J,fc=0 



(11) 



where \^° Bell ) = (|00) + 1 1 1 > ) / ^2 , S ab {t^ = -tr^N log 2 Xab{tJ] = 0, and (spin) 
negativity A/"[xas(t«)] = 1- Consider a density matrix xab and its partial transposition 
Xab f° r a two spin- 1/2 system AB. ya B is entangled if and only if Xab nas an y negative 



eigenvalues 



0Q- 



The negativity 



14| is a computable measure of entanglement defined 



by M[xab] = max{— 2 J]j ^i, 0}, where ^ is a negative eigenvalue of Xab- At r /> we wm 
show that 

I, , ,,,2 



X^(r/) = \jl £[R A 

z j,k=0 



S[Ri k ] 



1 + \u\ 



l + \u\ 2 
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1 + \u\ z 



i + M 2 



(12) 



with S AB (T f ) = -Plog 2 P-(l-P) log 2 (l-P) but P = (l-\u\ 2 )/2 (see FIG. 1 and 2), and 
-^1xab( t ./)] = l^| 2 ( see FIG. 3 and 4). Generalization to a system of iV spin-1/2 particles is 
straightforward: 

1 



XA 1 -A N [n i 



|vj>° 

l^GHZ 



)a 1 -a n (^ 



GHZ\ 



E K 



R 



,jk 
A N i 



j,k=0 



I 1 

XA v ~A N (l-f) = - ]T £{Rm\ 
j,k=0 



£[ R A N 



(13) 



Here, \*% HZ ) = (|0 • • • 0) + |1 • • • 1) )/y/2 Q. 

In order to measure the effects described by Eqs.(7) and (12), we introduce a static 
observer at each spacetime point along the "trajectory" of the particle(s). Each observer is 
assigned a local inertial frame defined by the following convenient choice of vierbein e£(x): 

1 „, . 1 



el(x) = l,ef(x) 



,e z 3 (x) = 1, 



(14) 



with all the other components being zero. Furthermore, we demand that the particle(s) be 
moving with four-velocity 



// / \ / ! ^ srnn £ sin -Q 
v>r{x) = (cosh 4 , — y===—, 0, smh £ cos v) 



(15) 



Vl + f 

or four-momentum q^(x) = mu ti {x). Here, tanh£ = v (= constant < 1), i.e., £ is the rapidity 
in the local inertial frame, and < •& < tt/2. In order for the particle(s) to move in this 
way, which is not a geodesic motion, we must apply an external force. The acceleration due 
to this external force is given by a M (x) = m a (x)Vam m (^): 

a» W = (Binh 2 5 sin* *| h, yT77, o, 



sinh 2 £sin 2 ^^-ln JT+f), (16) 
oz v 



where 

cosh£- + sinh£cos % In = ^ In y/l + f(t-z). (17) 

The inverse of the vierbein e^(x) in Eq.(14) is given by e°(x) = l,e^(x) = a/1 + /', e^(x) = 
\/l — /, e^(x) = 1. The vierbein transforms a tensor in a general coordinate system x^ into 
that in a local inertial frame x a . For instance, 

q a (x) = e^(x)q^(x) = (m cosh £, msinh£ sin i?, 0, msinh£ cos (18) 

and similarly, a a (x) = e^(x)a M (x) yields 

a a (x) = (sinh 2 £ sin 2 tf— In Jl + f,F sinh £ sin 0, - sinh 2 £ sin 2 — In a/TT/) . (19) 



A straightforward calculation shows that the nonzero components of the spin connection 
^(x) = e a x (x)V^(x) are = ^(x) = V^+fd In y/T+J/dt, u^x) = -u 3 xl (x) = 

VT+fd\nVT+f/dz, cu° y2 (x) = cu 2 y0 (x) = V^Jdlny/T^J/dt, and u 2 y3 (x) = -u 3 y2 (x) = 

Suppose at proper time r the particle(s) is at x^ . After an infinitesimal proper 
time dr, the particle(s) moves to a new local inertial frame at the new point x ffl = 
x v _|_ ?t f(| Ti q a (x) changes to q a (x') = q a (x) + 5q a (x) = A a b (x)q b (x), where the infinitesi- 
mal local Lorentz transformation A a b (x) = 5 a b + \ a b (x)dr, with \ a b (x) = — \a a (x)q b (x) — 
q a (x)a b (x)]/m + xV^) an d xV^) = — ^(^V^W- For our case, we have A° 1 (x) = 
A 1 (x) = sinh 2 £ cos sin A° 3 (:r) = X 3 (x) = - sinh 2 f sin 2 #G(t, 2; £, and 

A 1 3 (x) = — A 3 x (a;) = — cosh£ sinh£ sin z; £, •&); where G(t, z; £, i?) = (sinh£ cosW/dt + 



cosh£<9/<9z) In yl + /(£ — z) = —F(t,z;£,-d). Corresponding to A. a b (x) is the infinitesimal 
local Wigner rotation W a b (x) = 5 a fe + (f a b (x)dT, where <^°o( x ) = y?°j(^) — <P % o( x ) — an d 
<p l j(x) = X l j(x) + [Ao(x)/cj — k l \j (x)]/(\/k ■ k + m 2 + m). Its spin-1/2 representation is 
DW(W(x)) = a + i[ip23(x)a l + ip 31 (x)a 2 + ipi 2 (x)a 3 }dT/2, with the identity matrix a and 
the Pauli matrices {a 1 , a 2 , a 3 }. It follows that (p 1 ^) = —G(t, z; £, -d)H(k a \ £, ■&), with 

H(h a ;S,&) = (l- ^^^^^^tanhe) cosh £ sinh £ sin 0. (20) 
V yk-k + m 2 + m J 

Hence, for a finite proper time interval, 77 — Tj, we have 

/^(W^/,^)) = exp [-^(fcVi.T/,^)] , (21) 
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where 



1 

2 



\n^/l + f(t f - z f ) - \n^l + f(U- z t ) 
p(t f -Zf)-f{U-Zi)\H{h*W). (22) 

Consequently, \i/))a evolves to 

W) A = J N(k a )d 3 kJ2C(k%X)D^ 2 \w(x f ,xA)\A(x f ,x i )k a ,X') A , (23) 



A, A' 



and similarly 1^)^ to |*')ab = / N{k a )N(p b )d 3 kd 3 pJ2\,\>,a,*> C{k a , A; p b , a) x 
Dy x 2 \w(x f ,x l ))\A(x f ,x t )k a ,\') A ®D ( X 2) (W(x f ,x t ))\A(x f ,x t )p b ,a') B . We obtain Eqs.(7) 
and (12) by writing \ip') A and \^) A b respectively as density matrices, and tracing over the 
momentum degrees of freedom. This completes what we set out to do. 

In summary, we have shown that the spin entropy of a single massive spin-1/2 particle 
may change under the influence of a passing gravitational wave. Interestingly, this change 
has a dependence on the shape of the wave [see Eqs.(8) and (22)]. In other words, by 
determining the entropy change, one could in principle deduce /. To measure this change, 
one could prepare an identical ensemble of many particles in the state \ip) [Eq.(4)] and 
subject them to an external force that produces the acceleration in Eq.(16). The observers 
at each spacetime point then select a subensemble of particles to determine as accurately 
as possible its spin state. The variation of the spin entropy with proper time can then be 
determined. We may also consider the same experimental setup for two- or iV-particle state 
l^) [Eq.(lO) or its generalization, which gives Eq.(13)]. In this case, we can analyze the 
entanglement properties of the resulting states. Specifically, we have Af[xAB( T f)] = M 2 - 

We have to emphasize that the above effect, even though nonzero, is extremely tiny, 
especially in the light that the height or amplitude A of a gravitational wave may be of the 
order of 10~ 21 . Consequently, \u\ 2 would be extremely close to 1. So, in order to measure 
such a minute effect, we need to "amplify" or "concentrate" it. Our preliminary analysis 
of the 3- to 7-particle states shows that although a passing gravitational wave may have 
a greater effect on the 3-particle state compared to a 2-particle one, the 4-, 5-, 6-, and 
7-particle states are surprisingly robust. Thus, it seems, by considering iV-particle states 
(with N > 4) does not help. Here, we turn to another well-known phenomenon in quantum 



information science, entanglement swapping Q|. Briefly, we analyze the negativity of the 
resulting two-particle state 

E MA 2 = ^ BlB2 [(I AlM ® W^kJXxW 7 "/) ® xW T /))L (24) 

where \W Bell ) = (<r< ® a )^) (i = 0,1,2,3) and = tr[(/ AlA2 <g> I^W^kJ) X 
(x'aiBi( t /) ® x!42B 2 ( r /))] = 1/4 is the probability of obtaining outcome i from the Bell 
basis measurement. Particles B\ and B2 become maximally entangled, but ^aIa 2 yields 
A/"[H^ Aa ] = \u\ 4 . This is therefore an amplification of the decoherence effect due to a 
gravitational wave. We repeat the procedure with x'A j B j { T f) m Eq.(24) replaced by Br^ A 
to obtain , which has negativity AffH^^J = \u\ 8 . It is not difficult to see how one 

can achieve Af[E^ A2 ] = \u\ n , with n the number of particles. Hence, instead of a direct 
measurement on the spin states, we subject the particles to the above cycles of entanglement 
swapping, obtaining a smaller number of pairs of particles with negativities, which differ 
appreciably from 1. 

In conclusion, we have established the specific conditions under which the spin entropy or 
negativity of massive spin-1/2 particles may change due to a passing gravitational wave. This 
very small change may be amplified via the above entanglement swapping scheme, and may 
be measurable. It is therefore, a potentially viable means of gravitational wave detection. 
More generally, our results demonstrate the exciting possibility of detecting measurable 
effects due to spacetime curvature using ideas and tools developed in quantum information 
science. Effects inclu ding those due to our expanding universe will be discussed in a longer 
paper in preparation 16j|. 
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